arXiv: 1505.04046vl [astro-ph.SR] 15 May 2015 


Astronomy & Astrophysics manuscript no. oph acc © ESO 2015 

May 18, 2015 


X-Shooter study of accretion in p-Ophiucus: 
very iow-mass stars and brown dwarfs ^ 

C.F. Manara^’**, L. Testi^’^’"^, A. Natta^’^, and J. M. Alcala^ 

* Scientific Support Office, Directorate of Science and Robotic Exploration, European Space Research and Technology Centre 
(ESA/ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands 
e-mail: cmanara@cosmos. esa. int 

^ European Southern Observatory, Karl Schwarzschild Str. 2, 85748 Garching bei Miinchen, Germany 
^ INAF/Osservatorio Astrofisico of Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy 
Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching bei Miinchen, Germany 
^ School of Cosmic Physics, Dublin Institute for Advanced Studies, 31 Fitzwilliams Place, 2 Dublin, Ireland 
® INAF/Osservatorio Astronomico di Capodimonte, Salita Moiariello, 16 80131 Napoli, Italy 

Received March, 24th 2015; accepted May, 4th 2015 


ABSTRACT 

We present new VLT/X-Shooter optical and near-infrared spectra of a sample of 17 candidate young low-mass stars and brown dwarfs 
located in the p-Ophiucus cluster. We derived the spectral type and extinction for all the targets, and then we determined their physical 
parameters. All the objects but one have M*<0.6 Mq, and eight have mass below or close to the hydrogen-burning limit. Using the 
intensity of various permitted emission lines present in their spectra, we determined the accretion luminosity and mass accretion rates 
(Macc) for all the objects. When compared with previous works targeting the same sample, we find that, in general, these objects are 
not as strongly accreting as previously reported, and we suggest that the reason is our more accurate estimate of the photospheric 
parameters. We also compare our findings with recent works in other slightly older star-forming regions, such as Lupus, to investigate 
possible differences in the accretion properties, but we find that the accretion properties for our targets have the same dependence 
on the stellar and substellar parameters as in the other regions. This leads us to conclude that we do not find evidence for a different 
dependence of M^cc with M* when comparing low-mass stars and brown dwarfs. Moreover, we find a similar small (<1 dex) scatter 
in the Macc-44* relation as in some of our recent works in other star-forming regions, and no significant differences in Macc due to 
different ages or properties of the regions. The latter result suffers, however, from low statistics and sample selection biases in the 
current studies. The small scatter in the Macc-M* correlation confirms that mass accretion rate measurements in the literature based 
on uncertain photospheric parameters and single accretion indicators, such as the Ha width, can lead to a scatter that is unphysically 
large. Our studies show that only broadband spectroscopic surveys coupled with a detailed analysis of the photospheric and accretion 
properties allows us to properly study the evolution of disk accretion rates in star-forming regions. 

Key words. Stars: pre-main sequence - stars: formation - brown dwarfs - protoplanetary disks - accretion, accretion disks - open 
clusters and associations: individual: p-Ophiucus 


1. Introduction 


The evolution of protoplanetary disks surrounding forming stars 
has been a major subject of study in recent years. The interest 
in this topic is driven by the fact that disks are the birthplace 
of planets. A clear understanding of the physical mechanisms 
driving the formation, evolution, and dispersal of disks is thus 
needed to constrain planet formation theories, and to explain the 
observed properties of our own solar system and of the wealth 
of exoplanetary systems discovered so far. 

During their evolution, the central young stellar object 
(YSO) and the surrounding disk interact through multiple pro¬ 
cesses, such as photoevaporation, stellar winds, and accretion of 
matter. The latter is a result of viscous processes happening in 
the disk, which drive its secular evolution (e.g.,|Hartmann et^ 


1998|l, while photoevaporation and winds/outflows are thought 


to play a major role in the final disk dispersal (e.g., [Alexander et| 
al. 2014| and references therein). A good understanding of these 


* This work is based on observations made with ESO Telescopes at 
the Paranal Observatory under programme ID 085.C-0876. 

** ESA Research Eellow 


processes is needed to describe how disks evolve, and thus how 
planet formation takes place. Star-disk interaction processes can 
be studied observationally through the strong signatures they in¬ 
troduce in the spectra of YSOs. Accretion shocks give rise to 
continuum excess emission in the UV (e.g., Valenti et al. 1993t 


Gullbring et al.[|1998||2000l|Calvet & Gullbring||1998[|Calvet et| 


akf 2000|l and the prominent emission of permitted hnes 


the whole spectrum (e.g., Muzerolle et al.j 1998a|b|c| 2003 


[Natta et ah] |2004|l, while winds are traced by various forbid 


den emission lines (e.g., Hartigan et al. 1995| Rigliaco et al. 


2013[ Natta et al. 2014|l. In recent years, new instruments have 


provided the possibility of studying these processes simultane¬ 
ously in large samples of objects. In particular, the X-Shooter 
spectrograph mounted on the ESO/VLT telescope (jVernet et al. 


20111 provides us the possibility of studying YSOs with un¬ 


precedented levels of detail, as it simultaneously takes spectra 
at medium resolution (R~5000-20000) over the complete wave¬ 
length range from ~300 nm to ~2500 nm. We are thus able to 
study all the signatures of accretion and winds in the optical and 
near-infrared spectra of YSOs at the same time. 
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A significant effort has been made recently to study YSOs 
with X-Shooter. The accretion process has been studied for 
objects located in various nearby star-forming regions (cr-Ori, 
Rigliaco et al. 2012[ Lupus, Alcala et al. 2014[ Chamaeleon, 


Manara et al. 


2015 I and at various phases of disk evolution (e.g., 
transitional disks; Manara et al. 2014[ ). These works, among var¬ 
ious results, showed that the mass accretion rates (Macc) depen¬ 
dence on the stellar mass (M*) has a significant smaller scat¬ 
ter (~0.4 dex) than previously observed ( [Alcala et al.| 2014| l, 
and that there are transitional disks with comparable accretion 
rates as less evolved disks ( [Manara et al.j |2014[ ). These works 
focused on very low-mass, pre-main-sequence (PMS) stars and 
up to solar mass YSOs, with only a few objects below the hy¬ 
drogen burning limit (e.g., Rigliaco et al.[ 201 lb[ [Stelzer et alT 


2013|l. These objects were located in regions with age ~2-3 Myr, 


or more. However, to fully test models of disk evolution it is im¬ 
portant to derive accretion properties of YSOs for a large range 
of stellar masses and age. The work we present here aims to 
enlarge the sample of very low-mass stars (VLMS) and brown 
dwarfs (BDs) studied with X-Shooter, focusing on objects lo¬ 
cated in the younger (age~l Myr) p-Oph embedded cluster. 

The p-Oph cluster is located at a distance d ~ 125 pc 
(Lombardi et al. 2008[ Loinard et al. 2008 1 , and it is still highly 
embedded in the parental cloud, with values of Ay up to 50 mag 
in the densest core (Wilking & Lada [1983| l, and with a rather 
high value of total-to-selective extinction ratio Ry ~ 5.6 (e.g.. 


Kenyon et al.] 

1998t[Chapman et al.| [2009[[McClure et al.|[2010| 

Comeron et a 

. 2010J. The total population of the surrounding 


Ophiucus region is estimated to be of around 300 YSOs (Evans 
jet al.[ [2009] ), with ~200 YSOs located in p-Oph itself. Given the 
large quantity of YSOs present in this cluster and its proxim¬ 
ity, it has been extensively studied in the past at various wave¬ 
lengths from optical, near- and far-infrared, millimeter, radio, to 


X-ray (e.g.. 

Greene & Meyer| 1995 Luhman & Rieke 

1999 

Wilking et a 

l.| 20051 Alves de Oliveira et al.[ 2010[ 2012 

2013 

Bontemps et al. 2001 1 McClure et al. 2010| Gagne et al. 

2004|l 


Various works aimed to understand the accretion properties 
early ages have targeted p-Oph. In particular, |Natta| 
hereafter N06) have studied the vast majority of the 
Class II YSOs in this cluster, selected from the ISO sample of 
[Bontemps et al.|f2001[ ). They used near-infrared emission lines 
(PayS and Bry) to determine the accretion rates, and a statistical 
approach, based on the assumption that the whole population is 
located on the Hertzprung-Russel diagram (HRD) on a single 
isochrone at 0.5 Myr, to derive stellar parameters for most of 
the target^ This work showed that accretion rates for BDs in 
p-Oph are typically higher than for objects with the same mass 
in Taurus, and that various strongly accreting BDs are found in 
this region. However, this result strongly depends on the already 
mentioned assumptions, which can result in incorrect estimates 
of the stellar parameters, extinction, or accretion for single ob¬ 
jects, and only more precise analyses can shed light on this issue. 
With the set of X-Shooter spectra we present here, we are able to 
derive stellar and accretion parameters with significantly higher 
precision. With this in hand, we aim at addressing the still open 
questions comparing accretion properties for BDs and VLMS in 
this region and in other samples observed with X-Shooter. 


of YSOs at 


et al. (2006 


' Their results were obtained using an old distance estimate of 160 
pc, and have been corrected for the newly determined distance and re¬ 
calculated with the same method, but assuming a single isochrone at 1 
Myr for the whole population hy [Rigliaco et al.j ( |2011a[ (. Here we re¬ 
port in Tables [PTj all the values derived using the most recent distance 
estimated for the cluster of 125 pc. 


The paper is organized as follows. In Section]^ we describe 
the sample, observations, and data reduction; in Section we 
then discuss the method used to analyze the spectra and the de¬ 
rived properties of the objects, together with their main spectral 
features. In Section]^ we determine the accretion properties of 
the targets, and discuss the differences with previous estimates 
for the same targets; the comparison of accretion properties 
with results from other regions is then carried out in Section]^ 
Finally, we discuss our findings and conclusions in Section]^ 


2. Sample, observations, and data reduction 

We report here about new data collected with the ESO/VLT X- 
Shooter spectrograph ( [Vernet et al. 20111 during Pr.Id.085.C- 
0876 (PI Testi). This instrument acquires spectra of the tar¬ 
gets from A ~300 nm to A ~2500 nm simultaneously, split¬ 
ting the spectrum in three parts, usually referred to as the 
UVB {AA -300-550 nm), VIS {AA -550-1050 nm), and NIR 
(AA -1050-2500 nm) arms. 

to cover as 


The sample was selected from the one of N06 


many of the objects with detected accretion and estimated stel¬ 
lar mass M*<0.1 Mq as possible. It comprises 16 Class II YSOs 
and one Class I YSO that are analyzed here for the first time, 
one transitional disk (ISO-Ophl96) that has been analyzed in 
[Manara et al.|P014) l, and three Class III YSOs. We focus on the 
analysis of the 16 Class II YSOs and the Class I YSO. All these 
targets were included in the ISO sample compiled by Bontemps 
et al. ([200 l[l and most of them have previously been studied 


spectroscopically (e.g., [Natta et al.[ |2002| [2004] [20061 [Wilking] 
[et al.[[T999l[2005[ see also Table [A'l[ for literature data on these 
targets). Their ISO number, other names, and coordinates are re¬ 
ported in the first three columns of Table [T] 

All the objects have been observed in service mode using 
the 1.0x11" slit in the UVB arm and the 0.9x11" slits in the VIS 
and NIR arms, which lead to a nominal resolution R = T/A/l- 
4350, 7450, and 5300 in the three arms, respectively. Different 
exposure times have been adopted for each target depending on 
their estimated fluxes. The dates of the observations and expo¬ 
sure times are reported in the fourth and fifth columns of Table[T] 
The late spectral type (SpT) of most of the targets and the high 
extinction of the region explain the low SNR of most spectra, at 
least in the UVB arm. The only spectra with any detected sig¬ 
nal in the continuum in the UVB arm are those of ISO-Oph032 
and ISO-Ophl23, although with low SNR. Some objects have 
also very low SNR in the VIS spectra at 4 < 600-700 nm, or 
even at longer wavelengths. The NIR spectra of all targets have 
a good SNR. We report the SNR for the various spectra at diffe¬ 
rent wavelengths in the last three columns of Table [Tj The object 
ISO-Oph072 was reported by |McClure et al. ( 2010[ l to be mem¬ 
ber of a binary system with a separation of 3.62". The X-Shooter 
spectrum analyzed here is the one of the primary component of 
the system, while the secondary was not included in the slit. 

Data reduction has been carried out using the X-Shooter 


pipeline (Modigliani et al. 20101 version 1.3.7 and the same 


procedure as in Alcala et al.[(2014l. For the targets observed on 


2010-06-01, 2010-08-20, and 2010-08-28, the photometric stan¬ 
dard observed during the night of observation is not supported 
anymore in the pipeline. The flux calibration of the spectra has 
thus been obtained using the photometric standard star observed 
the following night, which was always photometric. The flux ca¬ 
libration of the spectra reduced with the pipeline has been com¬ 
pared with available 2MASS photometry and rescaled to this 
using synthetic photometry on the spectra to match the photo¬ 
metric flux. The correction factors are larger than 1 and usually 
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Table 1. Sample, observing log, and SNR of the spectra 


Object/other name 

RA(2000) 
h :m :s 

DEC(2000) 

Obs. date 
YY-MM-DD 

Exp. Time 
[s] 

SNR at A= 

700 855 1300 

ISO-Oph023/SKSl 

16:26:18.821 

-24:26:10.52 

2010-08-20 

4x750s 

0 

4 

61 

ISO-Oph030/GY5 

16:26:21.528 

-24:26:00.96 

2010-08-20 

4x750s 

15 

71 

66 

ISO-Oph032 / GY3 

16:26:21.899 

-24:44:39.76 

2010-06-01 

4x750s 

12 

29 

65 

ISO-Oph033/GYll 

16:26:22.269 

-24:24:07.06 

2010-08-26 

4x1800s 

0 

2 

29 

ISO-Oph037 / LFAM3 / GY21 

16:26:23.580 

-24:24:39.50 

2010-08-21 

4x750s 

0 

2 

36 

ISO-Oph072/WL18 

16:26:48.980 

-24:38:25.24 

2010-04-06 

4x480s 

3 

17 

51 

ISO-Oph087 

16:26:58.639 

-24:18:34.66 

2010-08-28 

4x750s 

0 

2 

29 

ISO-Oph094 

16:27:03.591 

-24:20:05.45 

2010-08-28 

4x750s 

0 

1 

10 

ISO-Ophl02/GY204 

16:27:06.596 

-24:41:48.84 

2010-07-23 

4x750s 

12 

41 

68 

ISO-Ophl 15 / WLl 1 / GY229 

16:27:12.131 

-24:34:49.14 

2010-07-30 

4x750s 

0 

1 

37 

ISO-Ophl 17 / WLY2-32b / GY235 

16:27:13.823 

-24:43:31.66 

2010-06-08 

4x480s 

3 

17 

92 

ISO-Ophl23 

16:27:17.590 

-24:05:13.70 

2010-08-27 

4x750s 

58 

74 

83 

ISO-Ophl 60 / B162737-241756 

16:27:37.422 

-24:17:54.87 

2010-08-27 

4x750s 

1 

15 

62 

ISO-Ophl64/GY310 

16:27:38.631 

-24:38:39.19 

2010-06-08 

4x750s 

6 

31 

45 

ISO-Ophl65/GY312 

16:27:38.945 

-24:40:20.67 

2010-06-08 

4x750s 

0 

1 

27 

ISO-Ophl76/GY350 

16:27:46.291 

-24:31:41.19 

2010-08-21 

4x750s 

1 

16 

89 

ISO-Ophl93 / B162812-241138 

16:28:12.720 

-24:11:35.60 

2010-08-21 

4x750s 

2 

28 

24 


Notes. Exposure times are the same in the three X-Shooter arms (UVB, VIS, and NIR). Wavelengths for the SNR calculation are reported in nm. 


lower than ~1.5, and always less than ~1.9. Indeed, these spec¬ 
tra have all been obtained with slits whose sizes are comparable 
with the seeing, so the slit losses are small. The only spectrum in 
which the conjunctions between spectra from different arms are 
not very good is ISO-Ophl02. In this object, the quality of the 
spectrum in the last ~100 nm of the VIS arm is not good, and 
this results in a bad matching of the VIS and NIR arms. Finally, 
telluric line removal has been performed using the IRAI ^task 
telluric following the procedure discussed in Alcala et al. ( 2014| l 
and using standard telluric spectra obtained close in time and air 
mass to the targets. 


3. Stellar and substellar properties 

Here we discuss the method adopted to derive the stellar and 
substellar properties of the targets, and also the results obtained 
and main features present in the spectra. 


3.1. Method 


Spectral classification for objects located in p-Oph is a difficult 
process, as these are usually highly-extincted and large excess 
because of disk emission, is also present at near-IR wavelengths 
(e.g., Luhman & Rieke 1999[ [Wilking et al. 2005 [ |Natta et 


[ah] |2606| |. Thanks to the very large wavelength coverage of our 
spectra, however, we can derive the spectral type (SpT) and ex¬ 
tinction (Ay) of the targets using multiple spectral features si¬ 
multaneously, thus reducing the degeneracies present in this ana¬ 
lysis. Given that the spectra have very low or zero SNR in the 
UVB arm, the analysis cannot be performed using the method 
described in [Manara et al.| ( |2013b[ ), which simultaneously deter¬ 
mines SpT, Ay, and accretion luminosity (Lacc)^ htting the UV 
and visible part of the spectra with a grid of models that includes 
the contribution of photospheric and accretion-induced emis¬ 
sion, and reddening. Similarly, the analysis methods used by 


^ IRAF is distributed by National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 


Table 2. Grid of photospheric templates 


Object 

SpT 

Teff [K] 

Ref 

TWA9A 

K5 

4350 

1 

TWA6 

K7 

4060 

1 

TWA25 

MO 

3850 

1 

TWA14 

M0.5 

3780 

1 

TWA13B 

Ml 

3705 

1 

synthetic spectrum 

M1.5 

3600 

2 

Szl22 

M2 

3560 

1 

synthetic spectrum 

M2.5 

3500 

2 

TWA7 

M3 

3415 

1 

TWA15A 

M3.5 

3340 

1 

Sz94 

M4 

3270 

1 

S0797 

M4.5 

3200 

1 

Par-Lup3-2 

M5 

3125 

1 

S0999 

M5.5 

3060 

1 

synthetic spectrum 

M6 

3000 

2 

Par-Lup3-1 

M6.5 

2935 

1 

synthetic spectrum 

M7 

2900 

2 

synthetic spectrum 

M7.5 

2800 

2 

synthetic spectrum 

M8 

2700 

2 

synthetic spectrum 

M8.5 

2550 

2 

TWA26 

M9 

2400 

1 

TWA29 

M9.5 

2330 

1 

synthetic spectrum 

LO 

2200 

2 


References. 1. [Manara et al.|j2013a| i; 2. [Allard et al.jpOTTj i 


Alcala et al. ( 2014| l or Rigliaco et al. (2012 1 cannot be adopted 
here, as all these need to perform a fit of the excess in the UV 
part of the spectrum to derive Lacc, and are well suited for the 
SpT determination only for objects with low extinction. 

The procedure we use here is the following: We collect a grid 
of photospheric templates, which cover the SpT from K5 to LO, 
with a typical step of 0.5 spectral subclasses throughout the en¬ 
tire M subclass. These are chosen primarily from the Class III 
spectra of jManara et al.| ( |2013a| ), which are observed spectra of 
nonaccreting PMS stars obtained with the same instrument. As 
this grid of templates is incomplete for SpT later than M6.5 and 
earlier than M9, and has no spectra with SpTs Ml.5, M2.5, and 
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M 6 , to complete our grid we make use of the synthetic spectra 
BT-Settl of |Allard et al.| ( |MTT] l, which are smoothed to match the 
resolution of the X-Shooter spectra. We adopt synthetic spectra 
with values of logg=3.5 (in cgs units), typical of young objects in 
p-Oph (e.g., Comeron et al. 2010| l, and effective temperatures 
(Tetf), corresponding to the given SpT according to the SpT-reff 
relation of |Luhman et al.| ( [2003[ ). The details of the adopted grid 
of templates is reported in Table The synthetic spectra repro¬ 
duce well, in general, the observed templates down to Tefr ~3050 
K, i.e., SpT~M5.5, while there are some discrepancies at lower 
temperatures. This, unfortunately, is also the region of our ob¬ 
served template library, which has a sparse sampling of spec¬ 
tral types. We describe the typical reddening toward the objects 
with the reddening law of Cardelli et al. (1989 1 using a value of 
Rv—5.6. 

For each target, we find the stellar parameters (SpT, Ay) with 
an automatic procedure that compares the observed spectrum to 
the various templates, which are artificially reddened (Ay ranges 
from 0 to 18 mag in step of 0.1 mag) and then normalized to 
the observed spectrum at A ~1025 nm. The comparison is done 
between the value of the mean flux across a window of ~4 nm in 
28 different points of the spectrum from A ~ 700 nm to A ~1725 
nm, which cover various molecular features particularly strong 
in the spectra of VLMS and BDs and sensitive to T^g, such as 
TiO, VO, and other molecular features (see |Manara et al.[|2013a| 
for a list) in the VIS arm, and H 2 O bands in the NIR arm. We 
limit our comparison to points within the aforementioned spec¬ 
tral range to avoid regions with large veiling due to accretion 
and low SNR spectra due to high extinction (A <700 nm), and 
regions where the contribution of the disk emission is substan¬ 
tial (/I >1700 nm). The best fit is derived by minimizinga^^ 


distribution, similar to the procedure by [Manara et al. 


on the two free parameters SpT and Ay. This is defined as 


(2013b I, 


2 

Alike 


=z 


fobsi^i^ fmodi.'^i) 
^obs(A[) 


( 1 ) 


where the index i denotes the various points included in the 
analysis, fobs is the flux of the observed spectrum, f^od that of the 
template, and o-gi,^ is the standard deviation on the flux of the ob¬ 
served spectrum. Uncertainties in the estimated parameters (SpT, 
Ay) are derived from the distribution, considered as a Ax^ 
distribution. Points where the SNR on the observed target is less 
than 5 are excluded from the Xgj^g calculation. To have enough 
SNR in the VIS arm, we use spectra smoothed using the boxcar 
smoothing procedure included in the IRAF splot package. This 
procedure convolves the spectrum with a rectangular box, whose 
width we set to 7 pixels, and results in a broadening of the nar¬ 
rower absorption and emission lines and features, but preserves 
their fluxes and equivalent widths. The results of the fit do not 
change using the smoothed spectra, but only the uncertainties on 
the results and the numbers of points included in the fit. At the 
same time, the different resolution of the templates does not in¬ 
fluence the result of the fit, as all the points where the Xg]^g is 
calculated are free of narrow absorption lines. We also include 
a constraint on the accepted fits, namely that the observed spec¬ 
trum should have a larger flux than the template at T >2115 nm, 
i.e., in the K-band, a part of the spectrum not used otherwise in 
the analysis. 

Various checks have been carried out on this procedure to 
determine its reliability, which we list here: i) Starting from our 
own grid of templates, we artificially add to their spectra vari¬ 
ous amounts of veiling due to accretion, of emission from the 


disk, and of reddening due to extinction. By applying the proce¬ 
dure just described in most cases we retrieve their correct stellar 
parameters. These stellar parameters only have big discrepan¬ 
cies in cases where the veiling due to accretion is substantial, 
i.e., the ratio of accretion to stellar luminosity (Lacc/T*) is larger 
than 0.5. This is expected to be the case for very few low-mass 
Class II YSOs (e.g., [Alcala et al.j [2014[ l. ii) We fit several X- 
Shooter spectra of Class II YSOs in Lupus by Alcala et al. ( [2014| l 
and retrieve results compatible with theirs, similarly to case i , as 
these objects are not strongly veiled. Hi) We run our procedure 
using only synthetic spectra for the later type objects or varying 
the resolution of the synthetic spectra. These different grids gen¬ 
erally lead to differences in the derived parameters well within 
their uncertainties, iv) We try to use spectra smoothed both in the 
VIS and NIR arms, while in general we smoothed only the VIS 
arm, and we obtain small differences, v) We try to use only non- 
smoothed spectra, with which we get similar results, but with 
larger uncertainties. 

The method just described allows us to derive SpT and Ay 
for all targets and lead to good best fits. We show in Appendix [B] 


(Figs. |B.1[ - B .61 the reddening corrected spectra with their 


best-fit template, also including additional Spitzer photometry 
from the literature. The values of 

from ~0.5 to ~3, with only two objects, ISO-Ophl02 and 
ISO-Ophl23, with>4 (see AppendixlX for further dis¬ 
cussion on these objects). We adopt a slightly ffiferent approach 
for the two objects with larger veiling due to accretion in our 
sample, ISO-Oph072 and ISO-Ophl23. For these two targets 
we exclude from the Xgy_g calculation all the points at A >1340 
nm, i.e., in the //-band. This choice is made as we expect that 
strongly accreting objects make a significant contribution to the 
emission at near-IR wavelengths from a disk, which is not in¬ 
cluded in our modeling. With this choice, we obtain the same 
result for ISO-Ophl23 as with the inclusion of all the points 
in the;^^^J^^ calculation, but with a lower value As we 

discuss in the following, the results for ISO-Ophl23 are highly 
uncertain and the derived stellar properties for this object should 
be used with caution. For ISO-Oph072, instead we derive a sig¬ 
nificantly later SpT (M3.5 instead of K7), but with a significantly 
smaller value ofxig^^ (2.2 instead of 4.7). We check the results 
also excluding the //-band for objects with strong infrared ex¬ 
cess due to an envelope, such as ISO-Oph037 and ISO-Ophl65 
(see Table A.l and Figs. B.l[[B.6[ l. For both of these targets, we 
derive Ay and SpT within the uncertainties of the normal pro¬ 
cedure including all the points. Flowever, as fewer points are in¬ 
cluded in the fit, the statistical significance of the results is lower. 
We decide then to use the values obtained using all points in the 
fit for these two objects, as well. 


3.2. Derived properties 


The stellar parameters derived from the best-fit results are re¬ 
ported in Tablej^with their Icr uncertainties. The first two panels 
of Fig.[2show the histograms of the SpT and Ay derived for our 
targets. The vast majority of the YSOs in our sample are of M- 
type SpT, with 7 out of 17 targets having SpT M6.5 or later. The 
derived values of Ay are in the range from ~0 mag to ~16 mag, 
and with a uniform distribution in this interval. We compare our 


findings with those in the literature in Sect. 3.3 


From the best-fit parameters (SpT, Ay) we derive Teff using 
the SpT-Teff relation of Luhman et al.[ ( [2003[ ). Stellar lumino¬ 
sity (L*) is derived from the ratio of the flux of the reddening 
corrected observed spectrum to the flux of the best-fit template 
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Table 3. Spectral types, extinction, and physical parameters of the p-Oph Class II YSOs 


Object 

SpT 

Teff 

Av 

L* 

R* 

M* 



[K] 

[mag] 

[Lol 

[Rg] 

[Mol 

ISO-Oph023 

M7 

2900)*’ 

Q 74-0.2 


0.79±0.20 

^•^^-0.01 

ISO-Oph030 

M7 

2900)™ 

4 5+^-2 

u.uu/_oooo 

1.03±0.04 

^•^^-0.01 

ISO-Oph032 

M6.5 

2935)™ 

0 6^^-^ 
^•^-0.1 

'-'•^D-J_0.000 

0.70±0.03 

^•^^-0.01 

ISO-Oph033 

M 8 

2700)*™ 

7 74-O.4 
-0.2 

0 005 

0.32±0.06 

^•^■^-0.01 

ISO-Oph037 

K7 

4060)™ 

16.1)°;^ 

A 41 C 4-0.000 

u.‘tiD_ooi7 

1.31±0.21 

^•^^-0.01 

ISO-Oph072* 

M3.5 

3340)™ 

8 . 2 )“] 

u.iuu_0 004 

0.97±0.04 

u.2)U-ooi 

ISO-Oph087 

M4.5 

3200)™ 

1 j.j_02 

u.iu:7_0 005 

1.08±0.05 

-0.01 

ISO-Oph094" 

M1.5 

3600)™ 

10 0 "^® ^ 
-0.1 


0.25±0.06 

0.40 

ISO-Ophl02 

M5 

3125+^^ 

7 74-0.1 
^•^-0.1 

0 047 

u.u^/_0 000 

0.74±0.02 

^•^‘^-0.01 

ISO-Ophll5 

M2 


15 

A 14^-1-0.006 

u.iH-^_0.035 

1.00±0.14 

^•*^^-0.07 

ISO-Ophll7 

M3.5 

3340)™ 

9 1+0.1 

0 991 +0.000 
'-'•^^^-0.000 

1.41 ±0.04 

^•J-J-O.Ol 

ISO-Ophl23* 

K7” 

4060** 

4.3" 

0.078" 

0.56" 

0.62" 

ISO-Ophl60 

M7.5 

2800)™ 

'-’•^- 0.2 


0.73±0.03 

^•^^-0.01 

ISO-Ophl64 

M 8 

2700)™ 

5 1+^2 


1.05±0.05 

0.05 

ISO-Ophl65" 

M2.5 

3500)[™ 

12 l"^®^ 
^^•^- 0.2 

0 034^^-®‘^ 
-0.001 

0.50±0.10 

A 4Q+0.08 
-0.08 

ISO-Ophl76 

M7.5 

2800)™ 


u.UD / -o ooo 

1.02±0.04 

0.06 

ISO-Ophl93 

M 6 

3000)™ 

7 4 + 0.2 
' • - 0.2 

^-0.003 

1.00±0.07 

^•^^- 0.01 


Notes. M, are derived using the evolutionary tracks of |Baraffe et al.| i n998| . * Object with strong veiling due to accretion. " These stellar parameters 
are very uncertain and should be used with caution. See Appendix[A|for a discussion on this target. ' Subluminous YSOs. M* reported with no 
uncertainties are for objects whose location on the HRD is at the edge of the tabulated evolutionary tracks. 



comparison of the position of the targets on the HRD with the 
evolutionary tracks by |Baraffe et al. ( T998) l. We also show the di¬ 
stribution of M* for our targets in the right-hand panel of Fig. [T] 


According to the parameters derived with the Baraffe et al. 
( 1998| l evolutionary tracks, six objects in our sample have a mass 
below the hydrogen-burning limit (ISO-Oph023, ISO-Oph032, 
ISO-Oph033, ISO-Ophl60, ISO-Ophl64, and ISO-Ophl76), 
and two more have a value of M* close to this limit (M*~0.1 
Mq, ISO-Oph030 and ISO-Ophl93). The majority of the YSOs 
have M*<0.6 Mq, and only one object has M*~lMo. 


Fig. 1. Histograms of the properties of the targets. The left-hand 
panel shows the SpTs, the middle panel shows the Ay, and the 
right-hand panels shows the M*. According to M* derived using 
Baraffe et al. ( |1998 1 evolutionary tracks, the sample comprises 
six BDs and two objects with mass just at the hydrogen-burning 
limit. 


at A ~1025 nm (K) and the known parameters of the template. 
In particular, when the best-fit template is a Class III YSO, we 
know its L* and distance ( [Manara et al.[ |2013a| ), so we can de¬ 
termine L* for our target simply considering the squared ratio 
of the distances and the ratio K of the two spectra, as in Manara 
et al. ( 2013b| l. On the other hand, when the best-fit template is 
a synthetic spectrum, L* is simply the total flux of the synthetic 
spectrum rescaled at the distance of p-Oph and multiplied by K, 
i.e., C*=47r-£/^_QpjjW-Ftot, template, which is a similar procedure as 
in jAlcala et al.| ( |2011| and references therein). Stellar radii (R*) 
are derived from the relation AttR^ctbT'^, while M* from 


The position of the targets on the HRD is shown in Fig. 
together with the Icr (red lines) and 3cr uncertainties (blue lines) 
on the estimated parameters. An additional <0.2 dex uncertainty 
on L* should be considered, as this is the error on the luminosity 
of the templates due to uncertainties in the distance and in the 
flux calibration of the spectra (e.g., Manara et al. 2013a| Alcala 
jet al.[|20T4| l. The HRD is shown with two different sets of evo¬ 
lutionary tracks, those by Baraffe et al.|( 199^ upper panel) and 


those by D’Antona & Mazzitelli ( 1994| bottom panel). The ma¬ 
jority of the targets (14 out of 17) is found to be located on the 
HRD at positions compatible with young (<5-10 Myr) ages ac¬ 
cording to both sets of models, and at ages <3 Myr according to 
|D’ Antona & Mazzitelli| ( |T994| ) evolutionary tracks. The latter set 
of models appears to reproduce better the positions on the HRD 
of the BDs in the sample, in particular, of ISO-Ophl76 and 
ISO-Ophl64. We however adopt the parameters derived using 
the evolutionary tracks from jBaraffe et al. ( 1998) 1 for consistency 
with previous studies of accretion in nearby star-forming regions 
(e.g., Alcala et al. |2014| l. 
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Fig. 2. Hertzsprung-Russell diagram for the p-Oph Class II 
YSOs analyzed here. The continuous lines show the isochrones, 
while the dashed lines show the low-mass pre-main sequence 
evolutionary tracks. The upper panel is obtained using the tracks 
by Baraffe et al. (1998 1 , while the bottom panel with tracks by 
D’Antona & Mazzitelli ( 1994|l. ISO numbers are reported for 


all the objects. Red lines represent the Icr error on the parame¬ 
ters, while cyan lines the her one. The objects ISO-Oph094 and 
ISO-Ophl65 are classified as being subluminous probably due 
to an edge-on disk for the former and to a still partly optically 
thick envelope in the latter, which is a Class I target. The stellar 
parameters and the position on the HRD of ISO-Ophl23 is very 
uncertain due to the very strong veiling due to accretion, which 
makes the spectral classification extremely difficult. 


For three objects (ISO-Ophl23, ISO-Ophl65, and 
ISO-Oph094) the position on the HRD is significantly closer 
to the main sequence, thus suggesting an older age than the 
rest of the targets. However, none of these targets are probably 
really old. The first object, ISO-Ophl23, is extremely veiled 
and this makes it difficult to detect the photospheric continuum. 
Therefore, the stellar parameters for this target are very uncer¬ 
tain. We discuss this case in detail in Appendix Regarding 
ISO-Ophl65, mid-infrared data have shown that this object is 
still embedded in its own parental envelope, and it is classified 
as a Class I YSO ( [McClure et al.| [201 Oj l. For this reason, the 
stellar parameters are uncertain. Ay is high, and an erroneously 


lower stellar luminosity estimate is possible. Finally, we suggest 
that ISO-Oph094 is a target whose circumstellar disk is 
edge-on, as its spectrum has an almost undetected continuum, 
also with a very low SNR in the NIR arm. We note that distance 
uncertainties for these objects are probably not leading to an 
erroneous location on the HRD, as a factor of ~2.5 larger 
distance would be needed to position ISO-Ophl65 closer to 
the bulk of the population, and this factor is even higher (~5) 
for ISO-Oph094. For the remainder of this work, we refer to 
ISO-Ophl65 and ISO-Oph094 as subluminous YSOs and use 
a different symbol in the plots to differentiate these from the 
other targets. Additional information and discussion on these 
objects is given in Appendix [A| 


3.3. Comparison with previous resuits 

Most of the targets analyzed here have been previously ob¬ 
served in near-infrared spectroscopic or photometric studies 


(e.g., Wilking et al.||199'9)|Luhman & Rieke]|1999||NaUa et ak 


2002 


et al. 


2004 


20061 or with optical spectroscopy (e.g., Wilking 


20051. The stellar parameters (SpT, Ay) derived in these 


works for our targets are reported in Table jA.lj and in Table |D.l 
forlN06l 

The comparison between the SpTs derived here and those 
obtained with spectroscopy in the literature (Table [ATT ] i is shown 
in Fig. where all the literature estimates are reported using 
circles for the most recent one and crosses for the older re¬ 
sults. All our SpT estimates agree within 1 or 2 spectral sub¬ 
classes with previous studies, with the only two exceptions be¬ 
ing ISO-Ophl23, which was classified as M3.5, while we find 
a very uncertain best fit with SpT K7, and ISO-Oph072, which 
we found to have a SpT M3.5, while it was previously classi¬ 
fied as K6.5 by [Wilking et al.| ( 2005| l. In both cases, the targets 
are strongly veiled, and thus their spectral type classification is 
subject to large uncertainties that can justify these differences. 
Finally, we derive for ISO-Ophl 17 a SpT of M3.5, which is in 
better agreement with the SpT M5 reported by McClure et al. 
( |2010| l than with the SpT K8 found by |Gatti et al.| ( [20()^ 

According to |N06| only four objects in our sample should 
have M* well above the hydrogen burning limit (ISO-Oph037, 
ISO-Oph072, ISO-Oph087, and ISO-Ophl 17). The result of 
our analysis confirms that all these YSOs are VLMS. At the 
same time, among the other 13 targets that were classified 
as BDs in N06 only 8 are confirmed as BDs by our ana¬ 
lysis, if we also consider ISO-Oph030 and ISO-Ophl93 as 
BDs. Three more targets that were classified as BDs are those 
with highly uncertain stellar parameters, as they are either 
subluminous (ISO-Oph094, ISO-Ophl65) or strongly veiled 
(ISO-Ophl23). Regarding the remaining two objects that were 
thought to be BDs, we classify ISO-Ophl02 as of M5 SpT and 
with a M* slightly above the hydrogen-burning limit (M*= 0.15 
Mq), and this estimate is also compatible with Wilking et al. 
( 2005| l. Finally, our derived M* for ISO-Ophl 15 is 0.5 Mq, 
while it was reported to be just above the hydrogen burning limit 
by N06 However, this object was reported to have a SpT MO by 
other authors ( [Gatti et ak] [2006 [ [McClure et al.[ |2010[ l, which 
confirms it is a low-mass YSO, not a BD. 

In Fig.|^we show the comparison between the values of Ay 
derived here with those derived by N06| for all the objects in our 
sample (cyan squares) and by other works for most of the targets 
(black circles). We see that most values agree within a typical 
difference up to ~1 mag. The largest differences are found for 
ISO-Ophl65, where both [N06 and McClure et al. (2010 1 report 
a value of Ay ~15-17 mag higher than our estimate Av=12.1 
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Fig. 3. Comparison of the SpT derived here for the targets 
(Sect. [3^ with those reported in the literature (see Table A.l for 
the references). Black circles report the most recent SpT estimate 
found in the literature, while black crosses are used for other SpT 
estimates available. The dotted lines show the +2 spectral sub¬ 
classes difference between the two values. The larger discrepan¬ 
cies are for ISO-Ophl23, classified here as K7 except with very 
uncertain parameters, and ISO-Oph072, which we classify as 
M3.5, and was reported to be K6.5 in the literature. 


Fig. 5. Comparison of the L* derived here for the targets 


(Sect. 1^ with those from 


N06 


discussed by Rigliaco et al. 


(20 


and corrected for distance as 


Ta|l. The dotted lines show 


difference between the L* values of a factor 2. The only object 
with a significantly larger difference is lSO-Oph033, but its L* 
has been revised by Comeron et al. ( 2010| l and the latter value, 
in agreement with ours, is reported with an empty symbol. One 
sigma uncertainties in our estimates are reported, and are some¬ 
times smaller than the symbol. 



ky [mag] - this work 


Fig. 4. Comparison of the Ay derived here for the targets 
(Sect. |3.2| l with those reported in the literature. The latter are 
shown with cyan squares in the case of values derived by |N06 
and in black circles when from other works (see Table |A.l| for 
the references). The dotted lines show the +1 mag difference be¬ 
tween the two values. One sigma uncertainties in our estimates 
are reported, and are sometimes smaller than the symbol. 


mag, and then for ISO-Oph094, reported to have Av=14.6 mag 
by N06 and found to have Ay^lO mag here. In one case the de¬ 
rived value differs from the estimates of lN06l but not with other 
studies (lSO-Oph033) or viceversa (ISO-Ophl 15), while the 
differences of Ay derived for lSO-Ophll7 are higher than 2 
mag with all previous studies, but in different directions. 

Finally, we compare the L* derived here with those deter¬ 
mined by |N06|and modified to account for the correct distance 


estimate (see Table |D3- This is a very interesting comparison, 
as L* is derived with two independent methods. Therefore, it is 
very important to note that the two methods lead to similar re¬ 
sults for all targets, with differences smaller than a factor ~2. 
The only exception is ISO-Oph033, but we note that the revised 
analysis by [Comeron et al. ( |2010| l leads to the same value of L* 
as the value we derived here. 

To summarize, we find, in most cases, values for Ay and L* 
to be compatible with those by |N06| and other literature studies, 
but that the derived M* differ from the literature estimate. The 
main reason for the different estimates of M* is the different Tefj 
determined here with a more robust method than those by |N06| 
and with better data than most previous spectroscopic studies. 


3.4. Spectral features 

We visually investigate the spectra with SNR>10 in the conti¬ 
nuum at ~700 nm for the presence of the Li I absorption line at 
4670.78 nm, a known proxy of youth in YSOs. We clearly detect 
this line in the spectrum of ISO-Oph032, while it is only ten¬ 
tatively detected in the spectra of ISO-Oph030, ISO-Ophl02, 
and ISO-Ophl23. In the latter, the SNR of the spectrum is very 
high, but the strong veiling makes all the absorption lines very 
weak. For all these targets, we report the presence of lithium for 
the first time, while the low SNR of the other spectra hinder the 
detection of this line in the rest of the sample. 

Various permitted emission lines are detected in the spectra 
of our targets, and these are used in the next section to derive 
the accretion luminosity of the targets. In few spectra with high 
SNR, we also detect forbidden emission lines. In particular, the 
[OI] A 6300 A line is detected in the spectra of ISO-Oph030, 
ISO-Oph032, ISO-Ophl02, and ISO-Ophl23. The analysis of 
these lines is out of the scope of this paper and will be discussed 
in a following work. 
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Fig. 6. Accretion luminosity derived from various emission lines 
luminosity for the p-Oph targets. Each subplot shows the value 
of log(Lacc/f'G) derived using the various indicators reported on 
the x-axis (CaK, Hd, Hy, Hy6, He^s^vam, Ha, He.i 667 nm, Pad, Pay, 
PayS, Bry) and in order of increasing wavelength. The red solid 
line is the average values obtained from the detected lines, while 
the dashed lines are the Icr standard deviation of this value. 
Black crosses are for upper limits out of the range of the plot, 
i.e., larger than —0.1. As upper limits depend on 

the rms noise of the spectra, they are significantly higher than 
measurements when the SNR is low, in particular, in the UVB 
arm. 


4. Accretion properties of p-Ophiucus young steiiar 
objects 

The accretion luminosity for our targets is derived using rela¬ 
tions between the luminosity of emission lines (Lune) and 
which have been calibrated by Alcala et al. (2014 1 . These repre¬ 
sent the only suitable way to determine Lacc from our spectra, as 
the UV-excess is not detected in the spectra because of the high 
extinction of the region. This is an indirect method that is less ac¬ 
curate than direct fitting of the UV-excess, but it has been shown 
that Lacc determined in this way are consistent with direct mea¬ 
surements of Lacc when multiple lines are used (e.g.,|Rigliaco et 
|ar] [M2l|M:alTet al.||M4| ). 

For this analysis, we then select 11 emission lines that are 
usually bright and whose Lune has a good correlation with Lacc, 
namely CaK, Hd, Hy, H/3, He^sgvnn,, Ha, He^ 667 nm, Pa5, Pay, 
Pa^, Bry. We determine the flux of the emission lines from the 
flux-calibrated and extinction-corrected spectra using an auto¬ 
matic Python procedure that determines the value of the conti¬ 
nuum and the extent of the line, which is then directly integrated 
with no additional modeling. All the line extent determination 
of the automatic procedures are checked by eye, and the derived 
flux of the lines are compatible with results obtained using diffe¬ 
rent methods, such as direct integration using the splot package 
under IRAK The error on the line flux is obtained propagating 
the Icr standard deviation on the continuum flux over the integra¬ 
tion window. For nondetected lines, we calculate the 3cr upper 
limits with the relationship 3 x FnoiseX^^, where F„oise is the rms 
flux-noise in the region of the line and Ad is the expected aver¬ 
age line width, assumed to be 0.2 nm. The fluxes and equivalent 
width (EW) of the lines, together with their respective errors, are 


reported in Table C.l 


All the objects analyzed here with detected Ha line have 
EWHff well above the threshold for accretors reported by White 
|& Basri| ( |2003| l, and thus are confirmed accretors. This line is 
not detected in the spectra of lSO-Oph087, lSO-Oph094, and 
ISO-Ophll5, which have SNR~0 on the continuum adjacent 
to the line. For the latter target, we are able to determine Lacc 
from the Pa/3 line and this leads to a value above the chro¬ 
mospheric noise level, i.e., the intensity of the chromospheric 
emission in a nonaccreting YSO, determined by [Manara et al. 
( 2013a| l for objects with this Tefj. The value of Lacc determined 
for ISO-Oph094 from the luminosity of the Pa/3 line seems com¬ 
patible with pure chromospheric emission. However, this object 
is probably observed edge-on, thus the measured line intensity is 
probably only a lower limit of the real line intensity. Finally, no 
emission lines are detected for ISO-Oph087, which could then 
possibly be a nonaccreting YSO. 

The luminosity of the lines is calculated as Lune = ■ 

/line, where d is the distance of p-Oph and /une the flux of the 
line. The value of Lacc is then determined as the average of the 
values of Lacc determined from the Lacc-Tune relations from the 
detected lines, while the error is the standard deviation of these 
values. Figure 1^ shows the values of Lacc determined with each 
line for all the targets. The accretion luminosity for the targets 
are reported in the second column of Table while the number 
of detected emission lines used to calculate Lacc is available in 
the last column of the same table. We estimate an upper limit on 
Tacc for ISO-Ophl87 from the upper limit on the Ha line and 
use this value in the analysis. 

Finally, Macc is determined from the values of Lacc just ob¬ 
tained and from the stellar parameters M* and R* determined as 
explained in the previous section. These paramete rs are related 


by the relation Macc= T25 ■ Laccf^*/(GM*) (e.g., Hartmann et 
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Table 4. Accretion luminosity and mass accretion rates of the 
p-Oph YSOs 


Object 

log i-acc 

[Lo] 

log Macc 
[Mo/yr] 

Detected lines 
[#] 

ISO-Oph023 

-3.51 

-9.85 

4 

ISO-Oph030 

-3.19 

-9.58 

7 

ISO-Oph032 

-4.09 

-10.49 

9 

ISO-Oph033 

-4.33 

-10.70 

1 

ISO-Oph037 

-1.20 

-8.49 

3 

ISO-Oph072 

-1.47 

-8.36 

8 

ISO-Oph087 

<-2.21 

<-8.92 

0 

ISO-Oph094 

-3.71 

-11.32 

1 

ISO-Ophl02 

-3.27 

-9.98 

8 

ISO-Ophll5 

-2.58 

-9.69 

1 

ISO-Ophll7 

-2.87 

-9.64 

3 

ISO-Ophl23 

-2.12 

-9.56 

11 

ISO-Ophl60 

-3.25 

-9.56 

4 

ISO-Ophl64 

-2.82 

-8.90 

8 

ISO-Ophl65 

-3.04 

-10.34 

2 

ISO-Ophl76 

-4.60 

-10.73 

1 

ISO-Ophl93 

-2.86 

-9.30 

3 


Notes. Lace is obtained as the mean of the values from the different emis¬ 
sion lines. Mace are derived using M* determined from the evolutionary 
tracks of |Baraffe et al.| l |1998^ . The last column reports the number of 
emission lines detected and used to derive Lacc- 



-3 - 2-10 1 

log(L*/Lg) 


Fig. 7. Comparison of accretion rate luminosity as a function 
of stellar luminosity for our sample with data from N06 for 
the same targets. Values from N06 are shown with red crosses, 
while our results with black circles. Empty black circles are used 
for subluminous targets and downward arrows for upper limits. 
Dashed lines are for different Lacc/C* ratios, going downward 
from 1, to 0.1, to 0.01, as labeled. The object most to the left- 
hand side of the plot in both samples is ISO-Oph033. 


|1998| ), which assumes accretion from a radius R,„ - 5/?*. 
The final value is reported for every object in the third column 
of Table m 


4. 1. Comparison with previous resuits 

The most comprehensive catalog of accretion rates in p-Oph was 
compiled by |N06| Here we compare the accretion properties we 
derived for our sample with the accretion properties compiled 
by N06 for the same objects. As already mentioned in Sect. |3.3| 
the results may differ in particular on the estimates of Tefj and 
M* due to the different methodologies. Indeed, Macc depends 
on various stellar parameters (M*, /?*, thus L* and Leff), which 
have been determined here with a more detailed analysis. As we 
reported in Sect. 1 |N06 did not have estimates for the SpT for 
most of the targets and had to assume that they were located on 
the HRD on a single isochrone after determining L* from the 
infrared colors. Our estimates of L* are consistent with those of 
N06, except that the Tgff estimated by N06| may differ substan¬ 
tially from ours in many cases, up to the extreme cases of the 
subluminous objects, which could not be classified as such by 


as the rest of the YSOs. The value of L* we determine here 
for ISO-Oph033 (bottom left point in both samples) is larger 


than that by N06 as discussed earlier, leading to a position of 
this target on the Lacc-T* plane closer to the bulk of the pop¬ 
ulation than before. The object with the lower Lacc/T* ratio is 
ISO-Ophl76, but its Lacc is still compatible with a line inten¬ 
sity being genuinely due to accretion, and not to chromospheric 
emission. Finally, the position of ISO-Ophl23 on the Lacc=0.1 
L* line is probably due to the already discussed uncertainties in 
its parameters. In reality, because of its strong emission, both 
in lines and continuum, this object is expected to fall closer to 
the Lacc=T* line in this diagram, as it was according to N06 
Appendix [A|provides further discussion on this target. 

_ The comparison of our measurements of Macc with those by 

|N06| is shown in Fig.[^ where we plot the logMacc vs logM* val¬ 
ues for the objects analyzed here, again using the black symbols 


N06 given their methodological assumption. This results in dif¬ 
fract M* and R* between the two works. Moreover, |N06| use 
the intensity of a single emission line to derive Lacc, while we 
measure this quantity using multiple, up to 11, emission lines. 
We are thus less affected by scatter in one single Lacc-Tiine rela¬ 
tion or by other processes that modify the intensity of a single 
line. 

We first show in Fig. the comparison of the Lacc-T* rela¬ 
tion for the objects included in our sample using our own results 
(black points) and those from |N0^ (red symbols). The two dis¬ 
tributions appear very similar, with the majority of the points 
located between the Lacc=0.1 L* and Lacc=0.01 L* lines. In gen¬ 
eral, most of the targets are located very close to the Lacc=0.01 
L* line, with slightly smaller dispersion than in the previous 
work reported here. This is also the case for the two sublumi¬ 
nous targets, which appear on this plot in the same position 


for our results and red symbols for those by N06 With respect 
to the results of |N06| the largest discrepancies are in the derived 
M*, which in many cases are larger in our work. This is a conse¬ 
quence of the later SpT determined for various objects with re¬ 
spect to those assumed by |N06| (see discussion in Sect. |3.3[ ). At 
the same time, the range of values of M^cc is similar for the two 
works, but our estimates lead to a smaller number of strongly 
accreting BDs than those of |N06| In particular, all the objects 
in our sample with reliable (sub)stellar and accretion parameters 
are found on the Macc-M* to follow the best-fit relation found 



more than an order of magnitude. This object appears to have 
a L* slightly higher than the rest of the objects analyzed here 
with similar SpT (see Fig. |^, and this could lead to a slightly 
higher Macc- The only targets that are located on the Macc-AL* 
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log(M*/Mg) 


Fig. 8. Mass accretion rate as a function of mass for the p-Oph 
sample. Values derived for the objects analyzed here are reported 
with black markers, values from N06 are shown with red sym¬ 
bols. Empty black circles are used for subluminous targets and 
downward arrows are used for upper limits. The continuous line 
represents the linear fit of this relation by |Alcala et al. (2014 1 for 
a sample of accreting objects in Lupus. The dashed lines repre¬ 
sent the 1 cr deviation from the fit. 



log(LJLg) 


Fig. 9. Accretion luminosity vs stellar luminosity for the whole 
p-Oph sample discussed here (black symbols), for the Lupus 
sample of Alcala et al. ( |2014[ ), reported with green symbols, and 

' ^( 2012) 1 and re¬ 


fer a sample in cr-Ori studied by Rigliaco et al. 


ported with yellow triangles. Upper limits are shown as down¬ 
ward arrows. Empty symbols are used for subluminous objects, 
and stars for strongly veiled objects. Dashed lines are for diffe¬ 
rent Lacc/T* ratios, going downward from 1, to 0.1, to 0.01, as 
labeled. The solid lines is the best fit for the GTO sample by 
INatta et al.| ( |20T4| . 


plane significantly below the bulk of the population are the two 
subluminous targets. 

We can thus conclude from this analysis that, in our limited 
sample of BDs and VLMS, we see a smaller number of heavily 


accreting targets than N06 This strongly suggests that BDs and 


VLMS in p-Oph accrete at a rate similar to those in other re¬ 
gions, as we discuss in the next section, contrary to the claim of 
N06 As the typical Lacc/T* ratio in our targets is similar to that 
determined by |N06| our suggestion is that the differences in the 
Macc-4^* relation is mostly because of the different (sub)stellar 
parameters determined here. Given that we do not have assump¬ 
tions about the positions of the targets on the HRD, and we have 
larger wavelength coverage that lead to a better constraint of Ay, 
SpT, and L*, we are confident that our results are more reli¬ 
able and that the distribution we observe in the Macc-4^*plane 
for our sample resembles the real distribution. A larger study on 
the whole sample of N06| with a data set similar to ours would 
further reinforce our finding. At present we do not attempt to 
fit our results as we have few objects sparsely sampled in the 


parameter space (L*, L^, 


c). 


5. Accretion in p-Oph compared with nearby 
star-forming regions 

As discussed in Sect. the main interests in the data set pre¬ 
sented here are to determine whether BDs in p-Oph have signif¬ 
icantly higher accretion rates than in other star-forming regions, 
and to study a very young region to compare the accretion prop¬ 
erties of its targets with those of targets located in older star¬ 
forming regi ons. W e have just discussed how our results differ 
from those of |N06| in the Macc-44* plane, leading to lower accre¬ 
tion rates for BDs than those reported by them. Here we compare 
our findings with those of Alcala et aLjl 20141 in the Lupus I and 
III clouds, and those by Rigliaco et al. ( |2012| l in the cr-Orionis 
(cr-Ori) cluster. These are selected as comparison sample as they 
include a significant number of VLMS and some (~5) BDs, and 


have been studied using spectra obtained with the same instru¬ 
ment. In addition, our accretion rates are determined using the 
Lacc-Tiine relation determined from these samples, hence the bi¬ 
ases in the Macc estimates are minimized than when comparing 
Mace estimates with other data sets. As these data have been col¬ 
lected as part of the Italian guaranteed time observation (GTO; 
Alcala et al. 201 l| l, in the following we refer to the whole com¬ 
parison sample as GTO sample. Our targets are also interesting 
as they enlarge the number of targets at low M*, in particular, 
there are six more BDs and two more objects with M* just at the 
hydrogen burning limit. 


5.1. Accretion luminosity and stellar luminosity relation 

We show in Eig.|^the position of our targets and those from the 
GTO sample (cf. legend in the plot and caption for symbols ex¬ 
planations) in the Lacc-T* plane. This is a particularly interesting 
relation as it is obtained independently from models of PMS evo¬ 
lution, thus it can be considered a “purely observational” plot. 
The plot also indicates the best-fit relation for the GTO sample 
by Natta et al. ( 2014[ ), in which the slope is 1.53+0.18. 

All the targets analyzed here are located in the same part of 
the Lacc-T* plane as those from the GTO sample. In particular, 
we do not have any extremely strong accretors (Lacc/T*>0.5), 
with the only possible exception being ISO-Ophl23, which is 
probably more intensely accreting than what is reported here. 
The number of objects between the L^cc- 0.01 L* and Lacc= 0.1 
L* lines decreases drastically at L*<0.1 Lq in the GTO sample, 
and the same appears to be true in our sample, as well. Also for 
this reason the best-fit relation has a slope significantly steeper 
than unity, implying that for higher luminosity targets the emis¬ 
sion due to accretion is relatively more important than in lower 
L* objects. This is also confirmed by the data collected here. 

The location of the two subluminous targets in our sample is 
similar to that of subluminous targets in the GTO sample (empty 
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log(M JMg) 

Fig. 10. Mass accretion rate as a function of mass for the whole 
p-Oph sample discussed here (black symbols) and for the Lupus 
sample of |Alcala et al. (20141, reported with green symbols. 
Upper limits are shown as downward arrows. Symbols are as 
in Fig. 1^ The continuous line represents the linear ht of this re¬ 
lation for the Lupus sample ( | Alcala et al. 2014| l, and the dashed 
lines represent the 1 cr deviation from the fit. 


symbols in the plot). Similar to what was found by Alcala et al. 
( |2014| l for the four subluminous targets reported here, we find 
that they appear in the upper boundary of the points in the Lacc- 
L* plane. 

Finally, our sample does not widen the spread on the Lacc- 
L* plane with respect to that found with the GTO sample. This 
suggests that the small spread in this relation found by Alcala et 


al. (20141 is not particular to the Lupus clouds, but is probably 


also common in other star-forming regions and is obtained when 
the data are analyzed with homogeneous methods. 

5.2. Accretion as a function of steiiar mass 

We present the values of logMacc vs logMacc for our and the 
GTO samples (cf caption and legend for explanation of the sym¬ 
bols) in Fig. [To] to gether with the best-fit relation derived by 
Alcala et al. (2014|l for the Lupus sample alone. The latter has a 


slope of 1.8+0.2 and is particularly interesting as it is based on a 
large sample, but has a very small scatter (cr=0.4 dex) from the 


best-fit relation compared to previous studies (e.g., Muzerolle et 

al. 

2003[ Mohanty et al. 2005| Natta et al. 2006[ Herczeg & 

Hi 

enbrandj |2008| [Rigliaco et a 

.[ [201 lal [Manara et al.[ |2012|l. 


Our sample does not substantially broaden the scatter of M^cc at 
any M*, whereas most targets appear to be consistent with the 
best fit of the relation and its dispersion. Only a cou¬ 

ple of BDs and the two objects with M* close to the hydrogen¬ 
burning limit in our sample seem to have a value of Macc at the 
upper edge of the dispersion from the best-fit relation reported 
by Alcala et al. ( |2014) l. However, the only one well above this re¬ 
lation is ISO-Ophl64, which, as discussed earlier, has a slightly 
higher L* than objects with similar SpT This suggests that BDs 
in p-Oph are not significantly stronger accretors than those in 
other star-forming regions, if they are strong accretors at all. As 
the majority of our targets follow the best-fit relation of Alcala 


et al. 


(2014 1 , we then suggest that the scaling of Mgcc with M* 


based on targets located in a region with higher density than in 
the Lupus clouds. 

As in the Lupus sample, the objects classified here as sub- 
luminous (empty symbols in the plot) are located in the M^cc- 
M* plane below the others. However, these objects are probably 
edge-on disks, thus they do not imply that the spread is larger 
as their stellar and accretion parameters cannot be constrained 
properly with our spectra. 

Assuming that the slope of the logM^cc-logM* relation is 
~ 1.8, we find that this is slightly shallower than that predicted by 
Padoan et al. ( 2005| l from models of Bondy-Hoyle accretion onto 
the forming YSO. This slope is instead the same as predicted by 


Dullemond et al.| (2006 1 assuming specific initial conditions, i.e., 
a small spread of angular momentum of cores at the first stages 
of star formation. The small spread of at any M* found by 


Alcala et al.|(2014}, and also present in this sample, would also 


suggest that the initial conditions at formation are similar for 
these objects. As one of the stronger predictions from the latter 
work is a very tight M^cc-Mdisk relation, we await future ALMA 
studies coupled with large samples of objects with reliably deter¬ 
mined stellar and accretion properties for a further confirmation 
of their theory. On the other hand, various studies suggested that 
the Macc-4^* relation is mostly a result of the evolution of disks 
under the effects of photoevaporation (e.g., [Clarke & Pringle| 
|2006[ [Ercolano et al.||2014) l. If this is the case, then our results 
seem to favor models of X-ray driven photoevaporation, which 
predict a slope of ~ 1.5-1.7 ( Ercolano et al.j 2014| l, with respect 
to those of EUV driven photoevaporation, predicting much shal¬ 
lower slopes of ~L3 5 ([Clarke & Pringlej |2006| l. As discussed 
also by Alcala et al.[ ( [2()14| l for the Lupus sample alone, there 
is in the data analyzed here no evidence of the double power- 
law behavior suggested by [Vorobyov & Basu| ( [2009| l, and the 
apparent bimodality of past data could be ascribed to mixing 
Macc determined with different methods and evolutionary mod¬ 
els. Einally, the upper envelope of our measurements has a slope 
similar to the bulk of the population and steeper than 1, contrary 
to the prediction of [Hartmann et al.| ( 2006[ l, which describes the 
Macc-Mi, relation as a consequence of layered disk accretion. To 
further constrain models using this relation, it is necessary to ap¬ 
ply methods similar to those used for the GTO sample or to the 
method described here when the UV-excess is not detectable, to 
analyze larger and complete samples of data to avoid poor statis¬ 
tics. A strong caveat to this discussion is, however, that the exact 
value of the slope is uncertain and depends on several parame¬ 
ters. In particular, the choice of the evolutionary model used to 
derive M* is crucial, as this can lead to different slopes, and even 
to different spread of Macc- Here we have elected to use one evo¬ 


is the same for VLMS and BDs. This relation is also indepen¬ 
dent of environmental conditions, as it is confirmed by our data 


lutionary model (Baraffe et al. 1998 1 to have results compatible 
with past analyses, but our results on the Macc-4f* relation can 
be model-dependent. 

Einally, we do not see in our objects, in comparison with the 
GTO sample, any clear evidence of dependence of Macc with the 
age of the regions. As just discussed, the objects analyzed here 
and those in the GTO are roughly located in the same position on 
the Macc-47* plane. This, however, can be the result of various 
effects. First of all, according to the similarity solution models 
by [Hartmann et al. (19981, one should expect ~0.5 dex difference 
in the values of Macc for objects with similar M* and with age of 
1 Myr or 3 Myr. This value is comparable with the uncertainties 
on our Macc estimates, and with the spread of values observed 
in a given star-forming region. Then, it is also possible that an 
age spread of few Myr is present in the ages of objects located in 
the same star-forming region. This would flatten any difference 
in Macc that might, in fact, be there. Therefore, we can only aim 
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to constrain the dependence of Macc with the age of the targets 
with signihcantly larger samples and more reliable estimates for 
the PMS stellar age. 


6. Conclusions 


We presented new observations of 17 VLMS and BDs with disks 
located in the p-Oph cluster. The spectra analyzed here were ob¬ 
tained with the VLT/X-Shooter instruments and were usually of 
good quality from ~750 nm to ~2500 nm, allowing us to deter¬ 
mine SpT and Ay for our targets by comparing the spectra in 
various molecular features with a set of photospheric templates. 
With this analysis, we determined the (sub)stellar parameters for 
all the targets. The total sample comprises six BDs, two objects 
with M* close to the hydrogen burning limit, eight VLMS and 
low-mass YSOs, and one solar-mass PMS star. 

We determined the accretion properties for the whole sample 
using the intensity of the emission lines present in our spectra 
and compared our findings with the most complete work avail¬ 
able in this held ( N06| ). Our results differ from those by N06 as 
we hnd a smaller number of strong accretors, in particular, in the 
BD regime. This suggests that BDs in this region are not accret¬ 
ing at a higher rate than in other star-forming regions, contrary 
to previous suggestions. 

A further constraint on this result comes from the compar¬ 
ison of our results with those in the Lupus clouds b y [Alcala et 
[ar] ( |2014] l and in the cr-Ori cluster by |Rigliaco et al.| ( |2012[ l, the 
GTO sample. Our objects follow the same Lacc-T* and Macc-A7* 
relations as those in the GTO sample. However, these hndings 
are based on very few objects. 

Finally, we do not hnd any evidence of differences of the 
accretion properties of targets located in the young (~1 Myr) 
p-Oph cluster or in the older (~3 Myr) regions targeted in the 
GTO survey. However, the differences expected are small and 
have various effects, such as observational uncertainties and age 
spread, and can make this analysis even more complicated. On 
the other hand, this work and other recent works that have used 
the capabilities of the VLT/X-S hooter spectrograph together 
with detailed analysis techniques (e.g., Rigliaco et al.[ 2012[ 
[Alcala et al.| |2014[ [Manara et al.j [2014[ ) have shown that the ob¬ 
servational uncertainties in the estimates of accretion rates and 
the spread of values can be signihcantly reduced with respect 
to the past. This suggests that it will be possible with complete 
samples in various star-forming regions to have a clear under¬ 
standing of the evolution of accretion with time and with the 
properties of the central PMS star. 
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Appendix A: Comments on individual objects 

Here we report the previous results available in the literature in 
Table |A.1[ These are discussed in Sect, 
our estimates. In the following, we discuss some peculiar objects 
in the sample. 

A.1. ISO-Oph094 

The location of this object on the HRD that we derive is that of 
a very old PMS object, much older than 100 Myr. However, the 
spectrum of this object has SNR~0 in the whole VIS arm and 


3.3 in comparison with 
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Table A.l. Data available in the literature for the targets included in this work 


Object/other name 

RA(2000) 
h :m :s 

DEC(2000) 

SpT 

Ay 

[mag] 

References 

McClure et al. i 

2010 

classitication 

ISO-Oph023 / SKSl / CRBR 2317.3-1925 

16:26:18.821 

-24:26:10.52 

M7.5,M5.5,M7 

10.0 

1,2,3 

Disk 

ISO-Oph030/GY5 

16:26:21.528 

-24:26:00.96 

M5.5,M6,M6 

2.8 

4,1,3 

Disk 

ISO-Oph032/GY3 

16:26:21.899 

-24:44:39.76 

M8,M7.5 

0.0 

4,3 


ISO-Oph033/GYll 

16:26:22.269 

-24:24:07.06 

M8.5,M6.5,M8.5 

8.0 

5,1,3 


ISO-Oph037 / LFAM3 / GY21 

16:26:23.580 

-24:24:39.50 

M0,K5 

16.1 

7,6 

Disk/FS 

ISO-Oph072/WL18 

16:26:48.980 

-24:38:25.24 

K6.5 

10.4 

4 

Disk 

ISO-Oph087 

16:26:58.639 

-24:18:34.66 





ISO-Oph094 

16:27:03.591 

-24:20:05.45 

M3? 

20.1 

8 


ISO-Ophl02 / GY204 

16:27:06.596 

-24:41:48.84 

M5.5,M6 

0.5 

4,3 

Disk 

ISO-Ophl 15 / WLl 1 / GY229 

16:27:12.131 

-24:34:49.14 

M0,M0 

11.5 

7,6 

Disk 

ISO-Ophl 17 / WLY2-32b / GY235 

16:27:13.823 

-24:43:31.66 

M5,K8 

7.5 

7,6 

Disk/FS 

ISO-Ophl23 

16:27:17.590 

-24:05:13.70 

M3.5 

4.4 

4 


ISO-Ophl60 / B162737-241756 

16:27:37.422 

-24:17:54.87 

M6 

6.0 

3 


ISO-Ophl64/GY310 

16:27:38.631 

-24:38:39.19 

M8.5,M7,M6 

5.7 

1,2,3 

Disk 

ISO-Ophl65/GY312 

16:27:38.945 

-24:40:20.67 

M2 

15.0 

7 

Envelope/1 

ISO-Ophl76/GY350 

16:27:46.291 

-24:31:41.19 

M6,M6 

7.0 

3,7 

Disk 

ISO-Ophl93/B162812-241138 

16:28:12.720 

-24:11:35.60 

M6 

7.5 

3 



Notes. Spectral types, extinction, and mass accretion rates determined from the following studie s: l.|Wilking et al.|(|T999[ ) ; 2.p..uhman & Rieke| 
(119991; 3.|Natta et all 120021; 4. Wilking et al.l(120051; 5.|Comer6n et all(120101; 6.|Gatti et al.|(|20061, V.lMcClure et al.|(|2010| and references 
therein), 8. [Alves de Oliveira et a .|( |20121 


also has low SNR in the NIR arm. The SpT we have derived 


is compatible with the highly uncertain estimates by Alves de 
Oliveira et al. (|2012|l, who suggest that this object has probably 


a SpT earlier than M3 and strongly extincted {Ay-20.1 mag). 
Also, their near-infrared spectra had very low SNR and their 
estimates are tentative. Even if we allow for higher values of 
Ay in our ht, we obtain the same best ht as that we report in 
this manuscript. As there is no other information available in 
the literature for these targets from mid-infrared data other than 
Bontemps et al. (2001|l, we are not able to further constrain its 


evolutionary class, which was reported to be Class II using ISO 
photometry. With the spectra available and the information in 
our hand, we cannot further constrain the properties of this ob¬ 
ject. Our classihcation for this object as subluminous is also re¬ 
lated to the tentative detection of forbidden emission lines in its 
spectrum. However, given the very low SNR of the continuum of 
this object, the subtraction of sky lines is not trivial and residuals 
are present even after careful extraction of the spectrum. 


A.2. ISO-Oph102 

The value of for the best fit for this object is ~12.4. As 

discussed in Sect. the spectrum of this target is of very poor 
quality in the last ~100 nm of the VIS arm. We have therefore ne¬ 
glected some points in this region to avoid biases due to the poor 
quality of the spectrum. Even if the value of is high, the 

best-fit template reproduces correctly the spectrum in the other 
features at shorter wavelength and in the NIR. Einally, our esti¬ 
mates are in good agreement with the literature. 


A. 3. ISO-Oph123 


As we mentioned many times, this object is very peculiar and 
the analysis of its spectrum does not lead to a satisfying esti¬ 
mate of its stellar properties. This object has been studied by 


Wilking et al. (20051 who reported a SpT of M3.5 with no par¬ 


ticular comments, only that there were strong emission lines in 


the spectrum. Scholz (20121 has also noted that this object is ex¬ 
tremely variable in the near infrared, with variations up to ~0.9 
mag in A-band. The flux calibration of the spectrum analyzed 
here is compatible with the 2MASS photometry in J- and H- 
bands, and it was only correct by a factor 1.12 to account for slit 
losses to match these bands. The magnitudes at optical wave¬ 
lengths reported in the literature are, instead, slightly lower than 
the observed flux of this spectrum. The procedure adopted here 
for the rest of the targets does not result in a low value 
for this target. The best-fit estimate is also probably uncertain 
as this object seems a very strong accretor, and we discussed 
that our procedure fails if the target has Lacc/T*^0.5, as could 
be the case for this object. As the excess due to accretion could 
significantly change the shape of the continuum, it is not pos¬ 
sible to determine the real reddening of the spectrum without 
considering all the components, such as excess due to accretion 
or reddening. We tried to fit this spectrum using the procedure 
by |Manara et al.| ( |2()13b| l, as the SNR in the UVB arm, although 
very low, is not zero. However, this automatic procedure is not 
able to reproduce the spectrum in the B aimer continuum region 
without further constraints on some parameters of the fit. We 
have then constrained the value of Ay to be between ~5 mag 
and ~1 mag, as this is the range of values for which the spread 
of Lacc derived from the various emission lines present in the 
spectrum is minimum. With this constraint and forcing L* for 
this object to be that of a 1.2 or 2 Myr old YSO, we obtain a 
best fit with SpT M5. The accretion rates that would result from 
the fit of the UV-excess with these parameters are one order of 
magnitude higher than those derived from the emission lines, 
however, and would result in a VLMS {M^ - 0.1 7 Mq) accreting 
at Macc~8xl0“^ Afo/yr. As we show in Eig. 


C.l 


many emission 

lines in this spectrum show strong absorption features in their 
profiles, and this could suggest that their real flux is larger. We 
do not expect the difference with the real flux to be larger than a 
factor ~2-3, which still is not compatible with the estimates from 
the UV-excess fit. The profiles of these lines suggest that this ob- 


13 
















































Manara et al.: X-Shooter study of accretion in p-Ophiucus 


ject might have a powerful jet, but more analyses are needed and 
are out of the scope of this work. 


A.4. ISO-Oph165 


As reported by McClure et al. ( 2010) l, this target is a Class I ob¬ 
ject. Therefore, its photospheric continuum is strongly extincted 
by the envelope that still surrounds the central star and the esti¬ 
mates of L* are biased by this. This is the reason why this object 
appears to be subluminous on the HRD. 
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Appendix B: Best fit of observed p-Oph spectra 
with tempiates 
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Fig. B.6. Same as Fig. 
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Appendix C: Permitted emission iine profiies in the 
X-Shooter spectra 
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Fig. C.l. Line profiles of permitted transition lines of ISO- 
Ophl23. 






























Table C.l. Extinction-corrected fluxes and equivalent widths of Balmer, Paschen, and Bracket lines 
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Table C.2. Extinction-corrected fluxes and equivalent widths of helium and calcium lines 
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Appendix D: Stellar and accretion parameters by |Natta et al.| ( |2006) corrected for the most recent distance 
estimate 


Table D.l. Properties of p-Ophiuci young stellar objects from |N06] using distance 
of 125 pc. 


# 

Object 

(ISO#) 

Class 

logC* 

[Le] 

logTeff 

[K] 

logM* 

[Mo] 

Line 

EW(line) 

[A] 

logZ^iine 

[^o] 

logZ^acc 

[^o] 

logMacc 

[Mo/yr] 

1 

ISO-OphOOl 

II 

0.16 

3.58 

0.02 

Pa/? 

<0.7 

<-4.36 

<-1.93 

<-9.01 

2 

ISO-Oph002 

II 

-0.21 

3.54 

-0.24 

Pa/? 

2.8 

-4.14 

-1.63 

-8.55 

3 

ISO-Oph003 

II 

-0.13 

3.55 

-0.19 

Pa/? 

3.9 

-3.92 

-1.33 

-8.28 

4 

ISO-Oph005 

III 

0.12 

3.58 

-0.01 

Paj3 

<0.5 

<-4.57 

<-2.21 

<-9.27 

5 

ISO-Oph006 

II 

0.07 

3.58 

-0.03 

Pajl 

19.0 

-3.03 

-0.12 

-7.18 

6 

ISO-Oph009 

II 

-1.13 

3.48 

-0.99 

Pajl 

<0.3 

<-5.84 

<-3.94 

<-10.46 

7 

ISO-OphOll 

III 

-0.29 

3.53 

-0.31 

Pa/? 

<0.8 

<-4.61 

<-2.26 

<-9.14 

8 

ISO-Oph012 

II 

-1.07 

3.48 

-0.95 

Pa/? 

<0.2 

<-6.05 

<-4.23 

<-10.75 

9 

ISO-Oph013 

II 

-0.26 

3.54 

-0.28 

Pap 

2.1 

-4.31 

-1.86 

-8.76 

10 

ISO-Oph014 

III 

0.01 

3.57 

-0.09 

Pap 

<0.5 

<-4.71 

<-2.40 

<-9.41 

11 

ISO-Oph016 

III 

0.49 

3.66 

0.20 

Pap 

-5.7 




12 

ISO-Oph017 

II 

0.54 

3.67 

0.22 

Pap 

<0.8 

<-4.08 

<-1.55 

<-8.81 

13 

ISO-Oph018 

III 

0.06 

3.58 

-0.04 

Pap 

<0.5 

<-4.67 

<-2.36 

<-9.41 

14 

ISO-Oph019 

II 

0.55 

3.67 

0.23 

Pap 

<0.9 

<-3.97 

<-1.40 

<-8.67 

15 

ISO-Oph020 

II 

-0.02 

3.56 

-0.11 

Pap 

1.2 

-4.30 

-1.85 

-8.85 

16 

ISO-Oph023 

II 

-1.85 

3.45 

-1.36 

Pap 

1.8 

-5.63 

-3.65 

-10.09 

17 

ISO-Oph024 

II 

0.26 

3.60 

0.07 

Pap 

8.9 

-3.21 

-0.36 

-7.47 

18 

ISO-Oph026 

II 

-0.48 

3.52 

-0.48 

Pap 

2.4 

-4.27 

-1.80 

-8.58 

19 

ISO-Oph028 

III 

0.36 

3.61 

0.12 

Pap 

<0.6 

<-4.33 

<-1.89 

<-9.03 

20 

ISO-Oph030 

II 

-1.42 

3.47 

-1.14 

Pap 

0.3 

-6.03 

-4.20 

-10.68 

21 

ISO-Oph032 

II 

-1.71 

3.45 

-1.29 

Pap 

0.4 

-6.13 

-4.34 

-10.79 

22 

ISO-Oph033 

II 

-3.26 

3.33 

nan 

Pap 

<0.7 

<-7.22 

<-5.82 

<nan 

23 

ISO-Oph035 

II 

-1.17 

3.48 

-1.01 

Pap 

-0.6 




24 

ISO-Oph036 

II 

0.70 

3.68 

0.31 

Pap 

0.7 

-3.92 

-1.34 

-8.64 

25 

ISO-Oph037 

II 

-0.79 

3.50 

-0.78 

Pap 

3.0 

-4.62 

-2.28 

-8.87 

26 

ISO-Oph038 

II 

0.00 

3.56 

-0.10 

Pap 

<0.3 

<-4.87 

<-2.63 

<-9.64 

27 

ISO-Oph039 

II 

0.99 

3.71 

0.37 

Pap 

<0.5 

<-3.87 

<-1.26 

<-8.54 

28 

ISO-Oph040 

II 

0.53 

3.66 

0.22 

Pap 

12.7 

-2.83 

0.15 

-7.11 

29 

ISO-Oph041 

II 

-0.61 

3.51 

-0.61 

Pap 

<0.3 

<-5.42 

<-3.36 

<-10.06 

30 

ISO-Oph043 

II 

0.01 

3.57 

-0.09 

Pap 

<1.0 

<-4.40 

<-1.99 

<-9.00 

31 

ISO-Oph044 

III 

-0.15 

3.55 

-0.21 

Pap 

<0.7 

<-4.69 

<-2.38 

<-9.32 

32 

ISO-Oph046 

II 

-0.11 

3.55 

-0.18 

Pap 

16.5 

-3.29 

-0.48 

-7.43 

33 

ISO-Oph047 

III 

-0.40 

3.53 

-0.41 

Bry 

<1.0 

<-5.15 

<-1.73 

<-8.55 

34 

ISO-Oph051 

II 

-0.67 

3.51 

-0.66 

Pap 

<0.4 

<-5.39 

<-3.33 

<-9.99 

35 

ISO-Oph052 

II 

-0.79 

3.50 

-0.77 

Pap 

5.9 

-4.28 

-1.83 

-8.42 

36 

ISO-Oph053 

II 

-0.81 

3.50 

-0.79 

Pap 

<0.3 

<-5.60 

<-3.62 

<-10.20 

37 

ISO-Oph055 

III 

-1.82 

3.45 

-1.34 

Pap 

<1.0 

<-5.90 

<-4.02 

<-10.46 

38 

ISO-Oph056 

II 

-0.56 

3.51 

-0.55 

Pap 

<0.5 

<-5.13 

<-2.98 

<-9.71 

39 

ISO-Oph058 

III 

0.58 

3.67 

0.25 

Pap 

<0.5 

<-4.25 

<-1.78 

<-9.05 

40 

ISO-Oph059a 

II 

-0.40 

3.53 

-0.41 

Bry 

<0.6 

<-5.37 

<-1.93 

<-8.75 

41 

ISO-Oph059b 

II 

-0.29 

3.53 

-0.31 

Pap 

<3.0 

<-4.05 

<-1.51 

<-8.39 

42 

ISO-Oph062 

II 

0.55 

3.67 

0.23 

Pap 

1.4 

-3.79 

-1.15 

-8.41 

43 

ISO-Oph063 

II 

-0.73 

3.50 

-0.73 

Pap 

<0.3 

<-5.53 

<-3.52 

<-10.13 

44 

ISO-Oph064 

III 

0.04 

3.57 

-0.04 

Pap 

<0.7 

<-4.51 

<-2.14 

<-9.20 

45 

ISO-Oph065 

I 

-1.35 

3.47 

-1.10 

Pap 

<10.0 

<-4.56 

<-2.20 

<-8.69 

46 

ISO-Oph066 

III 

-0.67 

3.51 

-0.66 

Pap 

<0.8 

<-5.01 

<-2.81 

<-9.47 

47 

ISO-Oph067 

II 

0.21 

3.59 

0.04 

Pap 

10.0 

-3.24 

-0.41 

-7.49 

48 

ISO-Oph068 

II 

0.30 

3.60 

0.09 

Pap 

0.6 

-4.32 

-1.87 

-8.99 

49 

ISO-Oph069 

III 

-0.95 

3.49 

-0.88 

Pap 

<0.4 

<-5.55 

<-3.55 

<-10.10 

50 

ISO-Oph070 

II 

-3.34 

3.32 

nan 

Bry 

<1.2 

<-6.80 

<-3.22 

<nan 

51 

ISO-Oph072 

II 

-0.86 

3.49 

-0.83 

Pap 

36.0 

-3.55 

-0.83 

-7.39 

52 

ISO-Oph073 

III 

0.57 

3.67 

0.24 

Pap 

<0.4 

<-4.32 

<-1.88 

<-9.15 

53 

ISO-Oph074 

III 

-0.25 

3.54 

-0.28 

Pap 

-1.5 




54 

ISO-Oph075 

II 

-0.97 

3.49 

-0.89 

Bry 

<3.9 

<-4.76 

<-1.39 

<-7.93 

55 

ISO-Oph076 

II 

-0.76 

3.50 

-0.75 

Bry 

<3.0 

<-4.70 

<-1.33 

<-7.93 
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# 

Object 

(ISO#) 

Class 

logC* 

[io] 

logTeff 

[K] 

logM* 

[Mo] 

Line 

EW(line) 

[A] 

logZ^iine 

[^o] 

logZ^acc 

[^o] 

logMacc 

[Mo/yr] 

56 

ISO-Oph078 

II 

0.20 

3.59 

0.03 

PajS 

<0.4 

<-4.62 

<-2.29 

<-9.37 

57 

ISO-Oph079 

II 

-1.52 

3.46 

-1.19 

PajS 

<3.0 

<... 

<... 

<... 

58 

ISO-Oph082 

III 

-1.99 

3.44 

-1.43 

Paj6 

<2.0 

<-5.71 

<-3.77 

<-10.19 

59 

ISO-Oph083 

II 

-0.04 

3.56 

-0.12 

PajS 

4.5 

-3.78 

-1.14 

-8.13 

60 

ISO-Oph084a 

II 

-1.27 

3.47 

-1.06 

PajS 

<5.0 

<-99.00 

<-99.00 

<-99.00 

61 

ISO-Oph084b 

II 

-1.27 

3.47 

-1.06 

Bry 

<0.2 

<-6.48 

<-2.93 

<-9.43 

62 

ISO-Oph086 

II 

-0.60 

3.51 

-0.59 

PajS 

<0.8 

<-5.00 

<-2.79 

<-9.50 

63 

ISO-Oph087 

II 

-1.10 

3.48 

-0.97 

Paj6 

2.4 

-4.90 

-2.67 

-9.18 

64 

ISO-Oph088a 

II 

0.38 

3.62 

0.13 

PajS 

10.5 

-3.05 

-0.15 

-7.30 

65 

ISO-Oph088b 

II 

0.41 

3.62 

0.14 

Paj6 

7.0 

-3.18 

-0.33 

-7.48 

66 

ISO-Oph089a 

II 

-0.87 

3.49 

-0.83 

Pa^e 

<0.5 

<-5.39 

<-3.33 

<-9.89 

67 

ISO-Oph089b 

II 

-0.86 

3.49 

-0.83 

Bry 

<0.7 

<-5.67 

<-2.21 

<-8.77 

68 

ISO-Oph091 

III 

0.64 

3.68 

0.28 

Pa^e 

-0.9 




69 

ISO-Oph092 

II 

0.91 

3.70 

0.35 

PajS 

8.4 

-2.75 

0.26 

-7.02 

70 

ISO-Oph093a 

II 

-0.51 

3.52 

-0.51 

Paj6 

<20.0 

<-99.00 

<-99.00 

<-99.00 

71 

ISO-Oph093b 

II 

-0.51 

3.52 

-0.51 

Bry 

<0.7 

<-5.15 

<-1.73 

<-8.49 

72 

ISO-Oph094 

II 

-1.95 

3.44 

-1.41 

PajS 

5.0 

-5.33 

-3.25 

-9.67 

73 

ISO-Oph095 

II 

0.11 

3.58 

-0.02 

PajS 

<0.6 

<-4.55 

<-2.18 

<-9.24 

74 

ISO-Oph096 

III 

-0.47 

3.52 

-0.46 

PajS 

-0.9 




75 

ISO-Oph097 

III 

-0.29 

3.53 

-0.31 

Paj6 

-1.1 




76 

ISO-Oph098 

II 

-0.06 

3.56 

-0.14 

PajS 

<0.7 

<-4.62 

<-2.28 

<-9.26 

77 

ISO-OphlOO 

III 

-1.76 

3.45 

-1.32 

Bry 

<0.7 

<-6.21 

<-2.69 

<-9.14 

78 

ISO-Ophl02 

II 

-1.43 

3.47 

-1.14 

PajS 

2.0 

-5.22 

-3.09 

-9.57 

79 

ISO-Ophl03 

II 

0.43 

3.63 

0.16 

PajS 

<0.6 

<... 

<... 

<... 

80 

ISO-Ophl04 

III 

-1.36 

3.47 

-1.11 

Bry 

<0.7 

<-5.87 

<-2.39 

<-8.87 

81 

ISO-Ophl05 

II 

0.38 

3.62 

0.13 

PajS 

1.7 

-3.81 

-1.18 

-8.33 

82 

ISO-Ophl06 

II 

-0.49 

3.52 

-0.48 

Paj6 

<1.0 

<-4.82 

<-2.55 

<-9.33 

83 

ISO-Ophl07 

II 

-0.58 

3.51 

-0.57 

PajS 

<1.3 

<... 

<... 

<... 

84 

ISO-Ophl08 

I 

4.14 

4.47 

1.14 

Bry 

1.3 

-1.68 

1.39 

-6.60 

85 

ISO-Ophl09 

III 

-0.38 

3.53 

-0.39 

Bry 

<0.5 

<-5.19 

<-1.77 

<-8.60 

86 

ISO-OphllO 

II 

0.50 

3.66 

0.20 

PajS 

<0.3 

<-4.41 

<-2.00 

<-9.24 

87 

ISO-Ophlll 

III 

-0.62 

3.51 

-0.61 

Paj6 

<1.0 

<-4.93 

<-2.70 

<-9.39 

88 

ISO-Ophll2 

II 

0.29 

3.60 

0.08 

PajS 

17.9 




89 

ISO-Ophll3a 

III 

-0.21 

3.54 

-0.24 

PajS 

-4.0 

-99.00 

-99.00 

-99.00 

90 

ISO-Ophll3b 

III 

-0.21 

3.54 

-0.24 

Bry 

-5.6 

-99.00 

-99.00 

-99.00 

91 

ISO-Ophll4 

III 

0.92 

3.71 

0.36 

PajS 

<8.0 

<... 

<... 

<... 

92 

ISO-Ophll5 

II 

-1.15 

3.48 

-1.00 

Paj6 

2.7 

-4.91 

-2.68 

-9.19 

93 

ISO-Ophll6 

II 

0.12 

3.58 

-0.02 

PajS 

<1.0 

<-4.32 

<-1.87 

<-8.93 

94 

ISO-Ophll7 

II 

-0.48 

3.52 

-0.47 

Paj6 

2.4 

-4.44 

-2.03 

-8.82 

95 

ISO-Ophll8 

II 

-0.38 

3.53 

-0.39 

Pa^e 

<20.0 

<... 

<... 

<... 

96 

ISO-Ophll9 

II 

-0.38 

3.53 

-0.39 

Bry 

<0.7 

<-5.05 

<-1.64 

<-8.47 

97 

ISO-Ophl20 

II 

-0.35 

3.53 

-0.36 

PajS 

<2.5 

<-4.33 

<-1.89 

<-8.74 

98 

ISO-Ophl21a 

II 

-0.17 

3.55 

-0.22 

PajS 

<0.6 

<-4.76 

<-2.48 

<-9.41 

99 

ISO-Ophl21b 

II 

0.36 

3.61 

0.12 

Paj6 

<0.5 

<-4.40 

<-1.99 

<-9.13 

100 

ISO-Ophl22 

II 

-1.20 

3.48 

-1.02 

PajS 

<10.0 

<... 

<... 

<... 

101 

ISO-Ophl23 

II 

-1.31 

3.47 

-1.08 

PajS 

21.9 

-4.10 

-1.57 

-8.06 

102 

ISO-Ophl24 

II 

-0.24 

3.54 

-0.27 

PajS 

<10.0 

<... 

<... 

<... 

103 

ISO-Ophl25 

III 

1.97 

3.99 

0.50 

PajS 

<20.0 

<... 

<... 

<... 

104 

ISO-Ophl26 

III 

-0.12 

3.55 

-0.19 

Paj6 

<0.9 

<-4.56 

<-2.20 

<-9.15 

105 

ISO-Ophl27 

II 

-0.19 

3.54 

-0.23 

Bry 

<0.0 

<... 

<... 

<... 

106 

ISO-Ophl28 

II 

0.27 

3.60 

0.08 

Paj6 

<0.4 

<-4.58 

<-2.23 

<-9.34 

107 

ISO-Ophl29 

II 

0.06 

3.58 

-0.04 

PajS 

<30.0 

<... 

<... 

<... 

108 

ISO-Ophl30 

III 

-0.07 

3.56 

-0.15 

PajS 

-0.8 




109 

ISO-Ophl32 

II 

0.57 

3.67 

0.24 

Pa^e 

<0.9 

<-4.01 

<-1.46 

<-8.73 

110 

ISO-Ophl33 

III 

0.23 

3.59 

0.06 

PajS 

<0.9 

<-4.28 

<-1.83 

<-8.93 

111 

ISO-Ophl34 

I 

2.17 

4.11 

0.56 

Paj6 

<12.0 

<... 

<... 

<... 

112 

ISO-Ophl35 

III 

0.31 

3.61 

0.10 

PajS 

-1.3 




113 

ISO-Ophl36 

III 

-1.28 

3.47 

-1.06 

PajS 

<0.9 

<... 

<... 

<... 

114 

ISO-Ophl38 

II 

-1.84 

3.45 

-1.36 

PajS 

<1.0 

<-5.95 

<-4.09 

<-10.53 

115 

ISO-Ophl39a 

II 

-1.01 

3.48 

-0.92 

Bry 

<2.8 

<-4.97 

<-1.57 

<-8.11 
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# 

Object 

(ISO#) 

Class 

logC* 

[io] 

logTeff 

[K] 

logM* 

[Mo] 

Line 

EW(line) 

[A] 

logZ^iine 

[^o] 

logZ^acc 

[^o] 

logMacc 

[Mo/yr] 

116 

ISO-Ophl39b 

II 

-1.01 

3.48 

-0.91 

PajS 

<6.0 

<-99.00 

<-99.00 

<-99.00 

117 

ISO-Ophl40 

II 

0.42 

3.63 

0.15 

PajS 

3.7 

-3.51 

-0.77 

-7.93 

118 

ISO-Ophl41 

I 

2.08 

4.07 

0.53 

Bry 

0.9 

-3.11 

0.10 

-7.52 

119 

ISO-Ophl42 

II 

-0.04 

3.56 

-0.12 

PajS 

<0.5 

<-4.72 

<-2.42 

<-9.41 

120 

ISO-Ophl43 

I 

0.69 

3.68 

0.30 

Bry 

1.8 

-3.74 

-0.47 

-7.77 

121 

ISO-Ophl44 

II 

-0.08 

3.56 

-0.15 

PajS 

<1.4 

<-4.32 

<-1.88 

<-8.85 

122 

ISO-Ophl45 

I 

0.06 

3.58 

-0.04 

Bry 

2.0 

-4.02 

-0.71 

-7.77 

123 

ISO-Ophl46 

III 

-0.41 

3.52 

-0.42 

Bry 

<0.9 

<-4.97 

<-1.57 

<-8.39 

124 

ISO-Ophl47 

II 

0.51 

3.66 

0.21 

PajS 

<0.5 

<-4.31 

<-1.86 

<-9.11 

125 

ISO-Ophl48 

III 

-0.36 

3.53 

-0.36 

Paf3 

<0.5 

<-5.00 

<-2.80 

<-9.64 

126 

ISO-Ophl49 

III 

-0.11 

3.55 

-0.18 


-0.9 




127 

ISO-Ophl51 

II 

-0.43 

3.52 

-0.43 

Paf3 

<0.5 

<-5.08 

<-2.91 

<-9.72 

128 

ISO-Ophl52 

III 

0.04 

3.57 

-0.04 

Paf3 

<2.0 

<-4.07 

<-1.53 

<-8.59 

129 

ISO-Ophl53 

III 

-1.55 

3.46 

-1.21 

Bry 

<0.0 

<... 

<... 

<... 

130 

ISO-Ophl54 

II 

0.07 

3.58 

-0.03 

Paj6 

<1.0 

<-4.38 

<-1.95 

<-9.01 

131 

ISO-Ophl55 

II 

0.57 

3.67 

0.24 

Paf3 

4.0 

-3.30 

-0.49 

-7.76 

132 

ISO-Ophl56 

III 

-0.77 

3.50 

-0.77 

Pafi 

<0.4 

<-5.43 

<-3.38 

<-9.98 

133 

ISO-Ophl57 

III 

-1.10 

3.48 

-0.97 

Bry 

<1.0 

<-5.76 

<-2.29 

<-8.80 

134 

ISO-Ophl58 

III 

-1.59 

3.46 

-1.23 

Pa/3 

<1.0 

<-5.74 

<-3.81 

<-10.27 

135 

ISO-Ophl59 

I 

0.16 

3.58 

0.01 

Bry 

<0.0 

<... 

<... 

<... 

136 

ISO-Ophl60 

II 

-1.73 

3.45 

-1.30 

Pa/3 

3.3 

-5.26 

-3.15 

-9.60 

137 

ISO-Ophl61 

II 

-0.01 

3.56 

-0.11 

Bry 

<0.8 

<-4.67 

<-1.30 

<-8.30 

138 

ISO-Ophl63 

II 

0.26 

3.60 

0.07 

Pa/3 

2.0 

-3.87 

-1.27 

-8.37 

139 

ISO-Ophl64 

II 

-1.47 

3.46 

-1.16 

Pa/3 

0.8 

-5.63 

-3.66 

-10.13 

140 

ISO-Ophl65 

II 

-1.34 

3.47 

-1.10 

Pa/3 

9.6 

-4.56 

-2.20 

-8.69 

141 

ISO-Ophl66 

II 

0.22 

3.59 

0.05 

Pa/3 

3.3 

-3.64 

-0.95 

-8.05 

142 

ISO-Ophl67 

I 

2.27 

4.14 

0.58 

Pa/3 

<4.0 

<-2.01 

<1.27 

<-6.43 

143 

ISO-Ophl68 

II 

0.11 

3.58 

-0.02 

Pa/3 

1.3 

-4.19 

-1.69 

-8.75 

144 

ISO-Ophl69 

III 

-0.26 

3.54 

-0.29 

Pa/3 

-0.7 




144 

ISO-Ophl69 

III 

-0.77 

3.50 

-0.76 

Pa/3 

<0.9 

<-5.06 

<-2.88 

<-9.48 

145 

ISO-Ophl69 

III 

-0.26 

3.54 

-0.29 

Pa/3 

-0.7 




145 

ISO-Ophl69 

III 

-0.77 

3.50 

-0.76 

Pa/3 

<0.9 

<-5.06 

<-2.88 

<-9.48 

146 

ISO-Ophl70 

II 

-3.34 

3.32 

nan 

Pa/3 

<3.0 

<-6.68 

<-5.09 

<nan 

147 

ISO-Ophl71 

II 

-0.78 

3.50 

-0.77 

Pa/3 

<4.0 

<... 

<... 

<... 

148 

ISO-Ophl72 

II 

-0.60 

3.51 

-0.59 

Pa/3 

<1.0 

<-4.91 

<-2.68 

<-9.38 

149 

ISO-Ophl73 

III 

-0.21 

3.54 

-0.25 

Pa/3 

<3.0 

<... 

<... 

<... 

150 

ISO-Ophl74 

III 

-0.96 

3.49 

-0.89 

Pa/3 

-0.7 




151 

ISO-Ophl75 

II 

-1.86 

3.45 

-1.36 

Pa/3 

15.7 

-4.76 

-2.47 

-8.91 

152 

ISO-Ophl76a 

II 

-1.22 

3.47 

-1.03 

Pa/3 

<0.5 

<-5.65 

<-3.68 

<-10.18 

153 

ISO-Ophl76b 

II 

-1.21 

3.48 

-1.02 

Bry 

<0.6 

<-6.02 

<-2.52 

<-9.02 

154 

ISO-Ophl77 

II 

-0.43 

3.52 

-0.43 

Pa/3 

<0.4 

<-5.13 

<-2.97 

<-9.78 

155 

ISO-Ophl78 

II 

-0.73 

3.50 

-0.73 

Pa/3 

2.0 

-4.72 

-2.41 

-9.03 

156 

ISO-Ophl79 

III 

-0.88 

3.49 

-0.84 

Bry 

<0.7 

<-5.68 

<-2.21 

<-8.77 

157 

ISO-Ophl80 

III 

0.31 

3.61 

0.10 

Pa/3 

-6.3 




158 

ISO-Ophl81 

III 

-1.30 

3.47 

-1.07 

Pa/3 

-5.5 




159 

ISO-Ophl82 

I 

0.22 

3.59 

0.05 

Bry 

5.4 

-3.41 

-0.17 

-7.26 

160 

ISO-Ophl83 

III 

-0.18 

3.54 

-0.23 

Pa/3 

<0.7 

<-4.72 

<-2.42 

<-9.34 

161 

ISO-Ophl84 

III 

0.22 

3.59 

0.05 

Pa/3 

-0.7 




162 

ISO-Ophl85 

II 

-1.08 

3.48 

-0.96 

Pa/3 

<0.2 

<-5.96 

<-4.11 

<-10.62 

163 

ISO-Ophl86 

III 

-1.13 

3.48 

-0.98 

Pa/3 

<1.5 

<-5.16 

<-3.02 

<-9.53 

164 

ISO-Ophl87 

II 

-0.56 

3.51 

-0.55 

Pa/3 

5.6 

-4.12 

-1.61 

-8.34 

165 

ISO-Ophl88 

III 

-0.17 

3.54 

-0.22 

Pa/3 

-0.5 




166 

ISO-Ophl89 

III 

-0.86 

3.49 

-0.83 

Pa/3 

<1.5 

<-4.94 

<-2.72 

<-9.28 

167 

ISO-Ophl90 

II 

-1.65 

3.46 

-1.26 

Pa/3 

<5.0 

<... 

<... 

<... 

168 

ISO-Ophl91 

III 

-0.03 

3.56 

-0.12 

Pa/3 

<0.5 

<-4.74 

<-2.45 

<-9.44 

169 

ISO-Ophl92 

III 

-1.06 

3.48 

-0.94 

Pa/3 

<2.0 

<-4.98 

<-2.77 

<-9.29 

170 

ISO-Ophl93 

II 

-1.22 

3.47 

-1.03 

Pa/3 

1.8 

-5.20 

-3.08 

-9.58 

171 

ISO-Ophl94 

II 

-0.73 

3.50 

-0.73 

Pa/3 

<0.5 

<-5.30 

<-3.21 

<-9.83 

172 

ISO-Ophl95 

II 

-0.10 

3.55 

-0.17 

Pa/3 

0.8 

-4.53 

-2.17 

-9.13 

173 

ISO-Ophl96 

II 

-0.67 

3.51 

-0.66 

Pa/3 

3.2 

-4.44 

-2.03 

-8.70 
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Table D.l. Continued. 


# 

Object 

(ISO#) 

Class 

logC* 

[Le] 

logTeff 

[K] 

logM* 

[Mo] 

Line 

EW(line) 

[A] 

logZyiijie 

[^o] 

[^o] 

logMacc 

[Mo/yr] 

174 

ISO-Ophl97 

III 

-1.11 

3.48 

-0.97 

Pafi 

<1.0 

<-5.29 

<-3.19 

<-9.71 

175 

ISO-Ophl98 

III 

0.73 

3.69 

0.31 

Pafi 

<0.4 

<-4.12 

<-1.60 

<-8.90 

176 

ISO-Ophl99 

II 

-0.13 

3.55 

-0.19 

Pafl 

1.7 

-4.23 

-1.75 

-8.70 


M* and age are derived using the evolutionary tracks of Baraffe et al. ( 1998|l and assuming a single isochrone at 1 Myr as explained 


in 


Rigliaco et al. ( 201 la| l. 
















